Pulde-Ferrell-Larkin-Ovchinnikov superfluidity in one-dimensional optical lattices 
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Spin-polarized attractive Fermi gases in one-dimensional (ID) optical lattices are expected to be 
remarkably good candidates for the observation of the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) 
phase. We model these systems with an attractive Hubbard model with population imbalance. By 
means of the density-matrix renormalization-group method we compute the pairing correlations as 
well as the static spin and charge structure factors in the whole range from weak to strong coupling. 
We demonstrate that pairing correlations exhibit quasi-long range order and oscillations at the wave 
number expected from FFLO theory. However, we also show by numerically computing the mixed 
spin-charge static structure factor that charge and spin degrees of freedom appear to be coupled 
already for small imbalance. We discuss the consequences of this coupling for the observation of the 
FFLO phase, as well as for the stabilization of the quasi-long range order into long-range order by 
coupling many identical ID systems, as in quasi-lD optical lattices. 

PACS numbers: 03.75.Ss, 03.75.Lm, 71.10.Pm 



I. INTRODUCTION 

Multi-component attractive fermionic systems with 
unequal masses, densities or chemical potentials have at- 
tracted continued interest for many decades in several 
fields of physics ranging from high-energji^ to condensed 



matter ' and, more recently, atomic physic 



,4.5.6 



The 



interplay between pairing and density imbalance of the 
different fermion species leads to a rich scenario, which 
includes the possibility of various exotic superconduct- 
ing states 7 -. In this context, the Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) phase 8 , has recently attracted a 
great deal of interest from both the experimental and 
the theoretical community-i^i 3 -^^. In the FFLO phase 
Cooper pairing occurs between a fermion with momen- 
tum k and spin f and a fermion with momentum —k + q 
(q ^ 0) and spin j. As a result, the superconducting or- 
der parameter becomes spatially dependent. Originally, 
the most favorable systems for the observation of the 
FFLO phase were predicted to be clean superconducting 
films in the presence of an in-plane (i.e. Zeeman) mag- 
netic field, above the so-called Clogston-Chandrasekhar 
limit£. Nevertheless, despite the fact that the original 
prediction dates back to more than thirty years ago, the 
FFLO phase has been very elusive to detect. 

The experimental realization of interacting trapped 
Fermi gases with population imbalance^ has renewed 
the hope of observing the FFLO, thus stimulating an in- 
tense theoretical activity 6 .^. So far, most of the theoreti- 
cal analysis has focused on 3D cold atomic systems. How- 
ever, as in the case of solid-state superconductors, the 
region of phase diagram where the FFLO phase has been 



found to be stable is quite small 2 -^. On the other hand, 
quasi-one dimensional or strongly anisotropic systems 
(such as coupled chains, heavy- fermion, organic, high-T c , 
and CeCoIn 5 superconductors) are believed to be good 
candidates for the realization of the FFLO phase d 11 ' 12 . 
Since the dimensionality of cold atomic systems can be 
easily tuned, and indeed cold atoms have already been 
successfully trapped in ID geometries^ 3 ., it seems natural 
to consider these low dimensional systems as the ideal 
candidates to observe non-homogeneous superconductiv- 
ity of the FFLO type. 

Many important results are available on the properties 
of spin-polarized ID Fermi systems with attractive inter- 
actions, which have been obtained by different methods 
and techniques. These include the Bethe-^4 nsatz solu- 
tions of certain exactly solvable models like the Hub- 
bard or Gaudin-Yang model s 14 ' 15 ! 16 ] 17 , as well as dif- 
ferent types of numerical approaches such as density- 
matrix renormalization-group (DMRG) 18 ' 19 : 20 and quan- 
tum Monte Carlo (QMC)S1, or field theoretical tech- 
niques like bosonizatio n 12 ' 18 . Very recently Orsoi 4 - 
and Hu et al*^ have studied the phase diagram of 
harmonically-trapped ID polarized Fermi gases by com- 
bining the exact solution of the Gaudin-Yang model with 
a local-density approximation. Mean-field theory has 
also been applied by Liu et al£^, although it is known 
that it has a number of limitations^ in ID, particu- 
larly as far as paring correlations are concerned. DMRG 
has been employed by Feiguin et al^ and Tezuka and 
Ucda 20 , and QMC by Batrouni et al. 21 to investigate the 
pairing correlations in the spin-polarized ground state 
of the attractive Hubbard model in the presence of har- 
monic trapping. Previously, Yangi£ used bosonization to 
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study the pairing correlations and the phase diagram of 
a single ID Fermi system as well as an array of weakly- 
coupled ID Fermi systems in the presence of a Zeeman 
field. 

Yang's analysis is valid only close to a continuous 
magnetic-field-driven transition from a uniform BCS 
phase and an FFLO phase, which he assumed to be- 
long to the commensurate-incommensurate universality 
class 2 ^. Another important assumption in Ref. Il2j is that 
charge and spin degrees of freedom are decoupled at low 
energies for small polarization. However, this scenario 
does not apply to the Hubbard model away from half- 
filling 1 -^^ nor to the Gaudin-Yang model^, which are 
the relevant models for current ID cold atomic systems. 
As we show in this work, charge and spin degrees of 
freedom are indeed coupled already for small polariza- 
tion, which leads to important differences as compared 
to the scenario described by Yang. We also numeri- 
cally demonstrate that the pairing correlation function 
exhibits prominent oscillations with a wave number equal 
(up to finite size corrections, see below) to the difference 
of Fermi wave numbers, <7fflo = I^ft — fcpx | , as pre- 
dicted by FFLO theory in lDi£ and in agreement with a 
number of Luttinger-type theorem a 27 i 28 . Thus, the finitc- 
wave-number oscillations in the pairing correlation func- 
tion can be regarded as due to the excess of n-\ — un- 
paired majority-spin fermions 12 . This is because in ID 
the Fermi wave number and the density are proportional 
to each other: k^n — Tm a . On the other hand, this rela- 
tionship is no longer linear in dimensionality higher than 
one, and in this case the oscillations in the order parame- 
ter are related to the center-of-mass momentum q of the 
Coopers pairs. These observations seem to indicate that 
in the ID case it is not entirely clear whether there is 
a strict close parallelism with higher dimensional FFLO, 
and in some respects the system can be also understood 
as a a coupled Bose-Fermi mixture of spin-singlet pairs 
(the bosons) and unpaired fermion o 15 ! 26 . 

The paper is organized as follows. Section Ull presents 
the Hamiltonian that we use to describe the system of 
physical interest, while Section HTT1 reports and discusses 
our main numerical results. Our main conclusions are 
briefly reported in Section Hvl 



II. THE MODEL 

We consider a two-component mixture with a total of 
N fermionic atoms loaded in a ID optical lattice of L sites 
(the lattice constant is taken to be unity). The fermions 
are assumed to interact via attractive on-site interactions, 
whose strength can be tuned e.g. by means of a Fesh- 
bach resonance. Sufficiently away from resonance(s), this 
system is modeled by the attractive Hubbard model: 

L-l L 

H = -t^2 {c\ a ct+ia + H.c.) - U n^nei , (1) 
a-,i=i e=i 



where t is the hopping parameter, c| CT (q ct ) is the creation 
(destruction) fermion operator in the ^-site (£ € [1,L]), 
a =|, l the pseudospin-1/2 index (in experiments this 
labels the two different atomic hypcrfine states of the 
mixture), U > is the strength of the on-site Hubbard 
attraction, fi£ a — cJ^Qo-, and N a — ^2 a ntcr- In order 
to simulate the effect of an external trapping potential, 
open-boundary conditions (OBC) breaking translational 
symmetry will be used (these are indeed the most suit- 
able conditions for the DMRG treatment of the above 
model^). Our calculations are performed in the canoni- 
cal ensemble, and the results apply only to lattices away 
from half-filling, that is, when N / L. In the calcula- 
tions, the spin polarization 6 = (JV-f —N^)/(N^ +N^) was 
varied by decreasing 2Vj while keeping constant the num- 
ber of "background" up-spin atoms iVj, from iVj_ = iV-j- 
(the unpolarized case, i.e. 6 — 0) all the way down to 
Ni = (the fully polarized case, i.e. 6=1). 

In the unpolarized case (6 = 0), all fermions pair into 
spin singlets due to the attractive on-site interaction. 
This yields a gap to all spin excitations and therefore 
spin-spin correlations decay exponentially with distance. 
Singlet superconducting and charge-density wave corre- 
lations exhibit a slower decay (of power-law type in the 
ground state of a thermodynamically large system), be- 
ing singlet superconducting correlations the ones that 
dominate at long distances in systems away from half- 
filing 2 ^. The aim of this work is to study the nature 
of superfluidity for < 6 < 1 as a function of the di- 
mensionless ratio U/t (in the fully polarized case, where 
6 = 1, H describes a system of N = iVj noninteracting 
fermions). The expectation values (...) of all operators 
below are understood to be taken over the ground state 

of n. 



III. NUMERICAL RESULTS AND DISCUSSION 

Due to the OBC (or, in general, to any external po- 
tential that breaks the Bloch translational invariancc 
of the lattice) the spin-resolved site occupation profiles, 
n&r = (nta), exhibit Friedel oscillations. This is illus- 
trated in Fig. [T] In the unpolarized 6 = case the Friedel 
oscillations in nj are in-phase with those in uj, giving rise 
to large-amplitude atomic-density wave a 31 ' 32 in the total 
site occupation + n^. As it is clear from the top 
panel of Fig. [1] in the general 6^0 case the total site 
occupation displays N± maxima associated with the for- 
mation of Ni spin-singlet pairs that are delocalized over 
the lattice. In the bottom panel of Fig. [TJ we show the lo- 
cal spin polarization — n^i (which could be measured 
through phase-sensitive optical imaging 5 ). For small 6 
(see e.g. the plot for JVf = 20 and N± = 16) the local 
spin polarization displays JVf — N± maxima correspond- 
ing to the number of fermions that are left unpaired. 
With increasing 6 though the spatial dependence of the 
local spin polarization becomes more complicated: the 
amplitude of the oscillations in the bulk becomes indeed 
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FIG. 1: Top panel: DMRG results for the total site occupa- 
tion nt\ + nil as a function of site position I for a system with 
AT-j = 20 fcrmions in L = 60 lattice sites and U/t = 5. The 
number of down-spin fermions is JVj_ = 20, 16, 12, 6 and 2 (the 
corresponding spin polarization is 5 = 0%, 11%, 25%, 54%, 
and 82%). Bottom panel: the local spin polarization — n^j, 
as a function of I for the same system parameters as in the 
top panel. The thin solid lines are just guides to the eye. 



smaller, thus making it hard to clearly identify iVj — iVj 
maxima. These, however, are not distinctive and unam- 
biguous signals of FFLO pairing. 

We thus proceed below to present a study of pairing 
correlations: the model described in Eq. |T]), in fact, can- 
not sustain any true long-range order-22 in ID, i.e. the 
ground-state expectation value of the pairing operator 
Kg — bi\Cc\ is zero. In the unpolarized case and for 
an extended system, the correlation function of the pair- 
ing operator Cm = (AJA^/) decays with a power law 
\t — l'\~ x l K p at large distances, where 1 < K p < 2 is 
an interaction-dependent Luttinger-liquid dimensionless 



parameter—. In the top panel of Fig. [5] we illustrate our 
DMRG results for the spin-polarization dependence of 
C«'=l/2 at U/t = 5, which measures real-space super- 
fluid correlations between the site £' = L/2 (the center 
of the trap) and all the other sites. For 5 = the power- 
law decay of the C«'=l/2 for \i — L/2\ ^ 1 is clearly 
visible. For finite 5, instead, the superfluid correlator is 
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FIG. 2: Top panel: DMRG results for the pairing correla- 
tor C«'=l/2 as a function of site position I for a system with 
iVj = 20 fermions in L = 60 lattice sites and U/t — 5. The 
number of down-spin fermions is TVj = 20, 16, 12, 6 and 2 (the 
corresponding spin polarization is S — 0%, 11%, 25%, 54%, 
and 82%). Bottom panel: Interaction contribution to the 
momentum-space pairing correlator AC(k m ) as a function of 
wave number k m = 7rm/(L + 1). The symbol coding is the 
same as in the top panel. Note the well-defined peaks at 
9fflo. In the inset we show AC(k m ) for the unpolarized sys- 
tem with Nf = JVj. = 20, which shows a narrow peak at k%. 
The thin solid lines are just guides to the eye. 



characterized by a distinctive oscillatory character— and 
a very simple nodal structure with exactly N-\ — iVj ze- 
roes. We have carefully checked that the long-distance 
decay of C«/ = £/2 is still power-law, signaling quasi-long 
range superfluid behavior also at finite S. 

A careful analysis of the oscillatory character of 
C«'=l/2 can be done by means of the Fourier transform 
of the pairing correlator 



C(k m ,k m ')= 2J Ctt' ( Pm(£)'Pm>(t') 



(2) 



where <p m (t) = [2/{L + 1)]~ 1/2 sin (k m £) [with k m = 
7rm/(L + 1), m = 1 . . . L\ are the eigenstates of the 
hopping term in Eq. |1}. The mode with zero wave 
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number is excluded from the allowed k m values due to 
the OBC. The lowest energy mode corresponds to k\. 
The diagonal part of the matrix C(k m , k m >) will be sim- 
ply denoted by C(k m ) = diag{C(fe m , k m >)} = C(k m , k m ). 
In the bottom panel of Fig. [2] we plot the difference 
AC(k m ) — C(k m )~C^(k m ) between C(k m ) and its value 
in the noninteracting gas (i.e. at U/t = 0), C 1 - ' (k m )^. 

At 5 = C(k m ) possesses a very narrow peak at k\ 
(see inset in the bottom panel of Fig. [2]), signaling quasi- 
long-range superfluid order of the conventional BCS type. 
For a finite S, instead, C(k m ) has a local minimum at k\ 
and a single well-defined peak appears at a wave number 
Qfflo = h + |fc F t - fcpj.|) where k Fa = nN a /(L + 1) 
are the spin-resolved Fermi wave numbers. The peak at 
9fflo in the Fourier transform of the pairing correlator, 
which is a direct consequence of the simple real-space 
nodal structure illustrated in the top panel of Fig. [2J is 
a clear-cut signal of FFLO pairing. 

The DMRG data shown in Fig. prefer only to a single 
value of U/t = 5. We now turn to illustrate the depen- 
dence of the superfluid correlation functions on U/t. In 
the top panel of Fig. [3] we illustrate the dependence of 
AC(k m ) on the interaction strength U/t for a fixed spin 
polarization S — 25%. On decreasing U/t the quasi-long- 
range FFLO order (i.e. the height of the peak at <?fflo)) 
which is emphatically strong for large U/t, survives all 
the way down to the weak coupling regime. This can be 
quantified better by analyzing the size of the anomaly T 
at k m — <7fflO) which is measured by the difference be- 
tween left and right (discrete) derivatives of C(k m ) eval- 
uated at (/FFLO i 

T = C(<7fflo + + C(<7fflo - k\) - 2C(g FF Lo) ■ (3) 

In the bottom panel of Fig. [3] we plot T as a function 
of U/t < 5. In this range T decreases in a smooth fash- 
ion to its noninteracting value (i.e. T = 0) as U/t is 
decreased to zero. In other words, for every finite 6, 
C(k m ) tends uniformly and smoothly to its noninteract- 
ing value C^°\k m ) as U/t is decreased towards zero. For 
sufficiently large values of U/t the FFLO phase can also 
be characterized by the peak visibility defined by 
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FIG. 3: Top panel: Interaction contribution to the 
momentum-space pairing correlator AC(fc m ) as a function of 
k m , for iV T = 20, N L = 12 (S = 25%), and L = 60. The inter- 
action strength U/t is decreased from 50 to 0.2. In the inset 
we show C(km) for U/t = 5, 1, and 0. Bottom panel: The 
FFLO anomaly T in the superfluid correlator [see Eq. <[3j ] 
as a function of U/t for iVf = 20 and iVj = 18 (squares, 

5 = 5%), Ni = 12 (circles, 8 = 25%), and N l = 6 (triangles, 

6 = 54%). Inset: peak visibility v [see Eq. Q] as a function 
of 2 < U/t < 5 for the same system parameters as in main 
body of the figure (y is about 30% at U/t = 50). The thin 
solid lines are just guides to the eye. 



are calculated from the following equations: 



C(<7fflo) - C(ki) 



C(q 



FFLO) 



C(h) 



(4) 



This quantity is plotted in an inset to the bottom panel 

of Fig. m 

Before concluding, we would like to illustrate the 
behavior of the density-density, spin-spin, and mixed 
density-spin static structure factors, S nn (k m ), S mm (k m ), 
and S nm (k m ). These are defined by the sum over all 
frequencies of the corresponding dynamic structure fac- 
tors^ that can be in principle measured through Bragg 
spectroscopy or Fourier sampling of time-of-flight im- 
ages—. In practice, S nn (k m ) , S mm (k m ), and S nm (k m ) 



S„ n (k m ) = diag{J 7 .T '.[(h e he>) - (fit) (ne>)]} 
S mm (k m ) = diag{T .T .[(mim t ) - (rn t )(m t )]} , 
S nm (k m ) = Ai&giT .T .[(him t ,) - (n t )(m t )]} 

(5) 

where fi£ = + fii^ and rhi = — n^i. In Figs. 0]- 
[5] we show the dependence of S nn (k m ), S mm (k m ), and 
S n m(k m ) on U/t for a slightly asymmetric system with 
7V T = 20 and N t = 18 (S ~ 5%). We remind the reader 
that in the unpolarized 8 = case S nn (k m ) has a peak at 
k m = 2fcp| = 2fcpx that signals real-space atomic-density 
wave a 31 ' 32 . In the spin-polarized case, this peak splits 
into two peaks at 2/cpt an d 2fcpx- This is clearly visible in 
Fig. 2] in the static density-density structure factor (top 
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FIG. 4: Top panel: the density-density structure factor 
Snn(k m ) as a function of k m for iVf = 20 and iVj = 18 and 
different values of U/t = 0.4, 1, 2, 3, 4, and 5 (from bottom to 
top). Bottom panel: the spin-spin structure factor S m m(k m ). 
The symbol coding is as in the top panel. 



panel), which presents a double-peak structure slightly 
below k m = 37r/4 (2fcp| 2ir/3 and 2fcpj w 37r/5 for 
the system parameters in this figure). This double-peak 
structure is not so visible in the magnetic structure fac- 
tor S mm (k m ), most likely because magnetic correlations 
near 2/cpt an d 2fcpx are still quite suppressed by the su- 
perfluid correlations, at least in the weakly polarized case 
(in the unpolarized case they are completely suppressed 
by the pairing gap). From Fig.[5]we note that S nm (k m ) is 
non-zero even at small k m , thus indicating that spin and 
charge degrees of freedom are coupled at long wavelength 
even for a small imbalance. 



Experimental signatures of the FFLO phase 

The most direct way to detect FFLO pairing would 
be to measure the pairing correlation function C«<. We 
would like to remark here that this correlation func- 
tion is, at least in principle, measurable via interfero- 
metric schemes^ in which two atomic wave packets are 
coherently extracted from the gas at different positions 



.46 



e 0.015 
Co 




FIG. 5: The mixed density-spin structure factor S nm (k m ) as 
a function of k m for JVj = 20 and N\ = 18 and different 
values of U/t = 0.4, 1, 2, 3, 4, and 5. The symbol coding is as 
in Fig. [4] Whereas the results at small U/t (see, for example 
filled squares corresponding to U/t = 0.4) are likely to be 
somewhat affected by finite-size effects, which prevent the full 
development of a pairing gap, the results at larger U clearly 
show a non- vanishing weight of Snm(k m ) at small k m , thus 
indicating that spin and charge degrees of freedom are coupled 
at long wavelength. 



and then are mixed by a matter-wave beam splitter. 
The atom counting statistics in the beam splitter output 
channels has been shown^ to reflect the spatial depen- 
dence of Cu> ■ 

The oscillations of the pairing correlations will also 
leave a detectable signature in the noise correlation o 37 ' 38 , 
G n (k, k') = {hk, t nk',i) - (hk,i){hk',i), where hk,a mea- 
sures the number of fermions with momentum k and spin 
(7 in a time-of-flight experiment. With increasing spin- 
polarization, in fact, the peak at k = —k' = (k, 0, 0) 39 i 40 
[here (1,0,0) is the direction along the axis of the ID 
system] will shift to a finite relative momentum (see e.g. 
the work by Yang in Ref. [l(] and the very recent DMRG 
calculation by Liischer et aZ4i). 

However, it is worth pointing out that the strength 
of noise signal in a strictly- ID system will be strongly 
affected by finite-size and temperature effects. This is 
because in ID the order is not long-range but quasi-long 
range, and therefore the slowest decay exhibited by cor- 
relations (like the pairing correlations) is a power law. 
Thus, in order to enhance the strength of the experi- 
mental signal for FFLO, it would be desirable to couple 
many ID systems, as in a tight 2D optical lattice (ar- 
rays of "atomic quantum wires'')^, so that the quasi- 
long range FFLO order can become true long-range or- 
der. The phase diagram of many coupled ID systems 
has been worked out in Ref. Qjl, where the author found 
that, at small polarization 6, true long-range 3D FFLO 
order will occur when the Luttinger liquid parameter for 
the charge excitations, K p , is larger than 3/2. In such 
a case the low temperature properties of the system are 
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dominated by hopping of pairs (i.e. Josephson coupling) 
rather than by single-particle hopping (the latter would 
turn the system into an anisotropic Fermi liquid, which 
could in turn become unstable to the FFLO state under 
appropriate conditions^) . 

However, the analysis of Ref. LL2j assumed the transi- 
tion from the unpolarized to the polarized case to be- 
long the the commensurate-incommensurate universality 
class. This assumption implicitly neglects the coupling 
between charge and spin degrees of freedom at low ener- 
gies, which is known to modify the behavior of physical 
observables at the transitio n s 2 ^ In this work this cou- 
pling has been demonstrated to exist also at long wave- 
lengths by an explicit numerical evaluation of the mixed 
static structure factor S nm (k m ) in a weakly polarized 
system (see Fig. [5J. Thus, the phase diagram depicted 
in Fig. 1 of Ref. [l2j seems not appropriate for coupled 
ID Hubbard (or Gaudin-Yangi£) models, and the transi- 
tion to long-range order will not take place in general for 
K p = 3/2 and may in general depend on the system pa- 
rameters (i.e. U/t and the lattice filling for the Hubbard 
model). 

To the best of our knowledge, a quantitative phase di- 
agram of coupled ID systems lacking spin-charge separa- 
tion has not yet been calculated. Furthermore, it is worth 
noticing that in cold atomic systems with short range in- 
teractions, another important factor must be taken into 
account, namely the relative strength of the pair hop- 
ping when compared to the single-particle hopping. The 
strength of the latter is given by t± (t± <C ep, where 
£f is the Fermi energy, for the analysis based on cou- 
pled Tomonaga-Luttinger liquids to hold), where t± is 
the hopping amplitude between two neighboring ID sys- 
tems. However, in absence of long-range interactions pair 
hopping can be only generated by (virtual) single-particle 
tunneling events, which at lowest order, yield a (Joseph- 
son) coupling strength of order ij_/A CT , where A CT is the 
spin gap. The phase diagram predicted in Ref. [TH is the 
result of a calculation which only compares the scaling 
dimensions of the pair hopping and single-particle hop- 
ping operators, and thus does not take into account the 
microscopic details of the coupling between ID systems. 



Thus, the stabilization of FFLO long-range order by a 
weak coupling between ID systems in the FFLO phase 
(the extreme anisotropic limit that could not be accessed 
by the authors of Ref. |43T) does not seem easily achievable. 
In turn, the most likely scenario for arbitrary polariza- 
tion is that the single-particle hopping will control the 
physics at low temperatures, and the system will behave 
as a spin-polarized normal Fermi liquid, which, in turn, 
could become unstable towards 3D FFLO ordering^ un- 
der appropriate conditions. 



IV. CONCLUSIONS 

In summary, we have shown how ultracold spin- 
polarized two-component Fermi gases confined in ID op- 
tical lattices are FFLO superfluids whose pairing corre- 
lation functions are characterized by a power-law decay 
and a simple nodal structure. However, we have also 
shown that charge and spin degrees of freedom appear to 
be coupled already for a small value of the spin polariza- 
tion. Finally, we have commented on the impact of this 
coupling on the detectability of true long-range FFLO 
order arising from Josephson coupling between ID sys- 
tems. 
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